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ABSTRACT: 

A  sum  of  normal  acceleration  and  pitch  rate  appears  to  be  the  best  vari- 
able to  use  to  control  aircraft  in  the  longitudinal  axis.    The  C*-Criterion 
specifies  that  the  time  response  of  this  quantity  called  C*  must  fall  in  a 
prescribed  envelope  for  all  speeds  and  altitudes.    It  is  equivalent  to  re- 
quiring the  control  system  to  hold  the  coefficients  of  a  certain  equation, 
which  describes  the  aircraft's  short-period  motion,  fixed.    This  is  done  by 
using  feedbacks  with  variable  gains.    The  gain-changing  mechanism  is  found 
using  gradient  techniques. 

The  system  was  shown  to  be  practical  by  an  analog  simulation.  It  was 
found  to  be  tolerant  of  instrument  noise,  elevator  hysteresis,  and  other  com- 
plications not  accounted  in  the  analytical  derivations. 

The  study  strongly  suggests  a  modification  of  the  C*-Criterion.     It  is 
proposed  that  the  coefficients  in  the  expression  of  the  quantity  C*  should 
be  representative  of  the  aircraft's  parameters  rather  than  have  the  values 
currently  in  use. 

This  task  was  supported  by:    Office  of  Naval  Research 

Navy  Department 


NPS57RJ9061A 
4    June    1969 


FOREWORD 

Three  Master's  theses,  currently  in  progress,  are  extensions  of  the 
studies  outlined  in  this  report  and  will  use  it  as  reference.    The  theses  are: 

"An  Investigation  of  Oscillations  in  an  Adaptive  Aircraft  Control  System  under 
Large  Input  Commands,"  by  Ens.  L.  S.  Wisler,  USN. 

"A  Study  of  an  Adaptive  Aircraft  Control  System  in  a  Self-Organizing  Con- 
figuration, "  by  Lt.  H.   M.  Richarde,  USN. 

"The  Extension  of  an  Adaptive  Aircraft  Control  Concept  to  Helicopters,  "  by 
Lt.   J.   M.   Hood,  USN. 

Most  of  the  work  reported  was  done  while  the  author  was  employed  by 
Honeywell,  Inc.    The  analytical  derivations  are  presented  but  only  the  con- 
clusions derived  from  supporting  studies  with  an  analog  computer  are  given. 
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Adaptive  Control  for  Aircraft 

The  usual  method  for  accommodating  an  automatic  flight  control  system  to 
the  wide  variations  in  dynamic  characteristics  of  the  airframe,  with  changes  of 
airspeed  and  altitude  over  its  flight  envelope,  is  to  change  control-loop  gains 
with  measurements  of  air-data.    The  dynamic  pressure  or  the  Mach  number,  as 
estimated  by  an  air-data  computer,  is  generally  used  as  the  parameter  with 
which  the  gains  are  scheduled.    The  electronic  art,  developing  at  a  bewilder- 
ing pace,  is  used  to  provide  reliable,  high-performance  control  hardware.    All 
flight  control  systems  currently  in  production,  with  the  exception  of  an  adap- 
tive system  for  the  F-lll  and  some  fixed-gain  controllers,  use  air-data 
scheduling.    This  includes  the  system  for  the  giant  C-5A  and  flight  controls 
proposed  for  the  new  cycle  of  commercial  transport  aircraft.     However,  the 
engineering  analysis  required  for  design  is  extensive  and  the  final  adjustments 
to  tailor  the  control  system  to  the  aircraft  during  flight  tests  are  expensive,  so 
there  has  been  the  hope  that  methods  of  control  would  be  found  which  did  not 
rely  on  scheduling  with  air-data  but  changed  parameters  by  measurements  of 
dynamic  performance.    These  would  ultimately  lead  to  much  simpler,  more  reli- 
able and  more  universally  applicable  systems. 

A  large  number  of  such  schemes  has  been  proposed  in  the  last  fifteen 
years.    Some  have  been  studied  in  flight  tests  and  only  one  has  gone  into  pro- 
duction.   A  brief  survey  of  the  older  ideas  is  given  by  Blakelock,  Reference  1. 
Boskovich  and  Kaufmann,  Reference  2,  describe  the  evolution  of  a  successful 
high-gain,  wide  bandwidth,  model-follower  system.    The  Bendix  system, 
described  by  J.  Bernadyn,  et  al.  ,  Reference  3,  was  flown.    The  North  American 


SIDAC  Controller,  developed  by  Shipley  and  his  associates,  References  4  and 
5,  has  been  extensively  analyzed  and  completely  designed  but  has  been  held 
short  of  flight  test  because  of  lack  of  funds.    These  devices,  while  being  very 
successful  in  most  respects,  have  debilitating  characteristics. 

Controllers  which  detect  limit  cycles,  either  intentionally,  as  in  the 
Honeywell  device  for  the  X-15,  or  unintentionally,  as  in  the  General  Electric 
control  system  for  the  F-lll,  are  upset  by  air  turbulence  unless  special  pre- 
cautions are  taken  which  then  reduce  the  effectiveness  and  increase  the  com- 
plexity of  the  system.    We  conjecture  that  the  Bendix  OLAC  was  not  carried 
further  because  it  probably  did  not  control  adequately  at  all  flight  conditions. 
The  SIDAC  will  probably  show  unwanted  oscillations  under  large  commands. 
It  is  a  high-gain  model-follower  and  is  complicated.    Since  reliability,  ease 
of  maintenance,  effect  of  component  failures  and  results  of  battle  damage  are 
vital  considerations  besides  performance  in  control  system  evaluation,  the 
relatively  complex  adaptive  schemes  have  been  losing  favor.    Indeed,  some 
people,  for  example  Horowitz  in  Reference  6,  argue  that  linear  control  should 
be  adequate. 

A  somewhat  different  point  of  view  has  been  taken  by  R.  C.  Hendrick  in 
developing  a  system,  described  in  Reference  7,  currently  undergoing  flight 
testing  at  Patuxent  River.     His  approach  is  to  use  only  the  minimum  amount  of 
information  from  the  dynamical  motion  that  is  needed  for  the  single  control  sys- 
tem to  be  applicable  to  a  whole  class  of  aircraft  without  additional  tailoring. 

Our  study  reverts  to  the  goal  of  providing  a  response  that  is  uniform  for 
all  flight  conditions  without  using  measurements  of  dynamic  pressure,  altitude 


or  angle-of -attack.    Information  obtained  from  inertial  instruments,  such  as 
accelerometers  or  rate  gyros,  or  from  measurements  of  control  surface  posi- 
tions and  command  inputs  is  admissible.    We  persist  in  this  area  despite  the 
general  disenchantment  with  adaptive  controls  not  only  because  it  appears  that 
a  really  practical  solution  for  aircraft  is  close  at  hand  but  also  because  the 
study  of  the  problem  yields  basic  understanding  of  aircraft  controls  and  the  dif- 
ficult question  of  requirements  for  suitable  handling  qualities. 

The  Approach 

Our  scheme  is  very  close  to  that  used  for  SIDAC,   (4,  5).    The  modifica- 
tion is  based  on  the  observation  that  the  equation  for  handling-qualities  C*- 
Criterion  is  very  similar  to  a  basic  short-period  equation  for  the  motion  of  the 
aircraft.    A  modest  feedback  and  feedforward  with  variable  gains  holds  the  co- 
efficients of  the  equation  fixed.    The  C* -requirement  may  be  met  by  choosing 
these  coefficients  to  be  the  same  as  demanded  by  the  Criterion  or  by  adding 
a  fixed  outer  loop.    The  mechanism  for  varying  the  gains  is  found  by  a  gradi- 
ent calculation  similar  to  the  SIDAC -ana  lysis .    We  achieve  the  advantage  of 
a  low-gain,  narrow- bandwidth  system  which  is  very  insensitive  to  instrument 
noise,  bending  modes,  accommodates  the  primary  control  system  and  satisfies 
the  requirements  directly  rather  than  using  a  model -following  technique. 

The  SIDAC-System  identifies  parameters  and  uses  this  in  turn  to  adjust 
gains.    Since  the  accuracy  of  identifying  the  several  coefficients  depends  on 
the  frequency  content  of  the  motion  and  the  particular  flight  condition,  it 
appears  that  a  system  which  calls  for  the  required  response  directly  should 
have  an  advantage.    This  has  also  been  argued  by  Hofmann  and  Best,  Reference 


8 ,  in  a  fairly  similar  approach  to  control  of  the  lateral-directional  axes. 

The  methods  of  synthesizing  the  gain  changing  mechanism  have  been 
developed  by  many  authors.    We  studied  a  Lyapumov-function  technique,  re- 
ported in  Reference  9.    This  was  later  rediscovered  by  Parks,  Reference  10. 
P.   M.  Lion,  Reference  11,  has  generalized  and  unified  the  several  methods. 

The  C*-Criterion 

The  most  difficult  problem  in  flight  control  design,  besides  considerations 
of  making  the  system  invulnerable  to  component  failures,  is  in  deciding  on  what 
flying  characteristics  will  be  acceptable  to  the  pilot.    The  requirement  for  longi- 
tudinal response  which  best  fits  into  an  analytical  formulation  is  that  promul- 
gated by  Tobie,  Elliot  and  Malcom,  Reference  12.    The  criterion  is  that  the 
time-response  trace  of  a  quantity  called  C*,  for  an  abrupt  force  applied  to  the 
stick,  must  fall  within  a  certain  envelope,  shown  in  Figure  1.     C*  is  the  sum 
of  the  normal  force  applied  to  the  pilot's  seat  plus  a  constant  multiple  of  the 
angular  velocity  in  pitch.    Thus,  it  is  a  linear  combination,  with  constant  posi- 
tive coefficients,  of  normal  acceleration  at  the  aircraft's  center-of-gravity ,  the 
pitch  velocity  and  the  pitch  acceleration,  assuming  the  pilot's  station  is  ahead 
of  the  eg.    This  blended  response  variable  has  been  used  for  flight  controls  be- 
fore.    It  is  a  reasonable  quantity  to  consider  because  at  flight  conditions  with 
high  dynamic  pressures,  at  which  the  aircraft  is  very  responsive,  the  pilot 
likes  to  feel  the  normal  acceleration  with  stick  force.     He  bases  his  assess- 
ment of  the  aircraft's  performance  on  it.     However,  in  low-q  flight  conditions, 
for  example  in  landing,  pitch  rate  is  a  more  important  clue  to  the  aircraft's  con- 
trollability.   Indeed,  for  a  helicopter  at  hover,  pitch  attitude  is  the  variable 


which  must  be  controlled.    The  sum  of  the  quantities  exhibits  this  change  in 
emphasis.    The  steady-state  response  of  normal  acceleration  has  the  factor  of 
the  aircraft's  velocity  when  compared  to  that  for  the  pitch  rate  and  hence  nor- 
mal acceleration  will  dominate  at  high  speeds  and  pitch  rate  will  be  important 
at  low  speeds.    A  more  comprehensive  analysis  is  given  in  the  report  by 
Makers,  Reference  13. 

A  Conjecture 

The  study  of  our  system  by  analog  simulation  showed  that  it  was  neces- 
sary to  hold  the  coefficients  of  the  short-period  equation  fixed  as  values 
representative  of  the  ranges  of  the  aircraft's  coefficients  over  the  flight  enve- 
lope and  not  to  achieve  the  C*-values  directly.    Analysis  presented  later  also 
shows  this.     It  is  necessary  because  of  higher  order  effects  that  are  ignored 
in  the  basic  derivation  but  are  present  in  realistic  situations.    The  C*-criterion 
is  met  by  using  an  outer-loop  feedback  with  fixed  coefficients. 

This  naturally  prompts  a  review  of  the  C*-requirement.    For  the  F-4  air- 
craft, which  was  used  as  the  vehicle  for  the  study,  the  pilot  sits  only  twelve 
feet  ahead  of  the  center-of-gravity.    This  makes  the  effect  of  pitch  accelera- 
tion on  the  C*-quantity  very  weak.    Since  most  of  the  studies  of  C*  were  done 
on  fixed-base  simulators,  it  appears  possible  that  the  angular-acceleration 
clue  is  not  properly  accounted  in  the  criterion.    Some  studies,  for  example, 
done  by  Bihrle,  Reference  14,  indicate  that  it  is  an  important  clue.    Thus,  we 
conjecture  that  a  more  natural  criterion  for  longitudinal  handling-qualities  is 
that  the  coefficients  of  the  C*-expression  be  taken  in  the  range  of  the  coeffi- 
cients of  the  corresponding  short-period  equation. 


Not  only  does  this  make  our  adaptive  scheme  more  natural  but  it  appears 
that  the  envelope  of  allowable  responses,  Figure  1,  may  be  made  much  more 
restrictive.    Then  the  Criterion  becomes  more  definitive  and  easier  to  apply. 
Further,  it  can  be  directly  extended  to  include  helicopters  by  adding  a  term 
for  pitch  attitude.    The  corresponding  short- period  equation  for  the  helicopter 
in  hover  has  just  such  a  term  and  this  quantity  is  an  important  handling  cue  in 
this  condition.    Thus,  this  modification  of  C*  is  intuitively  appealing. 

Equations  for  Short-Period  Motion 

We  begin  with  the  equations  of  perturbations  from  straight  and  level  flight, 
written  as 

(i) 


in  terms  of  angular  rate  in  pitch  £,  angle-of-attack  of  ,  and  normal  accelera- 
tion at  the  center-of-gravity  ft .    The  quantity  Oq  is  the  value  of  the  aircraft's 
unperturbed  velocity  and  £    represents  the  elevator  deflection.    The  coeffi- 
cients  f^t  ,  £^    and  x*  are  constants,  representative  of  the  flight  condition, 
and  have  appropriate  dimensions.    All  angles  are  measured  in  radians  and 
linear  quantities  are  measured  in  feet.    Data  for  the  F-4  aircraft  is  given  in 
Table  1. 


Angle-of-attack  is  difficult  to  measure,  and  its  interpretation  is  compli- 
cated by  turbulence  in  the  air,  so,  following  Shipley  (4),  we  eliminate  it  from 
the  equations  and  find 


The  coefficients  are 

P.  =  -  MV2.  -  "V*  (3) 

and  are  listed  in  Table  1.    They  are  negative  for  all  flight  conditions.    The 
first  equation  in  Set  2  is  the  fundamental  relation  in  this  study. 
The  C*-quantity  is 

Csn+{&     +     oc£  (4) 

in  which  j£     is  the  distance  of  the  pilot  forward  of  the  eg  and  Uc  is  a  number 
called  the  cross -over  velocity,  usually  taken  around  four  hundred  feet  per  sec- 
ond.    If  we  require  C*  to  be  precisely  a  multiple  of  the  command  input  C^   , 


10 


C  S5      -     te   CC    j  (5) 


Equation  4  may  be  written  as 


•• 


°t   a       X  v.      -t  r  <5) 


This  is  exactly  the  form  of  the  fundamental  equation  in  (2). 

The  Basic  Derivation 

The  control  configuration  is  diagrammed  in  Figure  2.    To  avoid  complicat- 
ing the  discussion,  for  the  present  neglect  the  actuator,  actuator  model,  the 
fixed  feedback  and  the  dynamics  of  the  primary    control  system.    Thus,assume 
for  now  that 


(7) 


We  are  looking  for  a  control  function 

f    -    H  e    +    fl>  vi  +   Vt  C  (8) 


in  which  the  gains    fi       C_    (*L    vary  so  that  for  each  flight 


condition 
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the  over -all  system  behaves  as  if  the  aircraft  equation 


were  replaced  by  one  with  given  fixed  coefficients   pL  ,    By.  ,    ft- 


(9) 


i* 


e-^©-Pv,v>-p»c    =0.  (io) 

Hence,  we  define  an  error  signal,  on  which  the  generation  of  the  variable 
gains  is  based,  by  the  amount  with  which  Equation  10  fails,  namely 

£  x.    b  —  p^e   -  p**vi  -  foe.  (id 

Subtracting  (9)  from  (11)  and  incorporating  the  control  function  yields 


(12) 


Therefore,  to  hold  the  error  zero  the  gains  must  be  adjusted  to  the  values 
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Values  of  these  with 

p\=-3    ,     Pv,»-.oa.  t    p,  =  -20 

are  listed  in  Table  1. 

To  construct  the  mechanism  for  calculating  the  gain,  we  follow  Shipley  in 
considering  the  gradient  of  the  sum  of  squares  of  the  coefficients  of  Equation  (12). 
This  method  is  now  a  standard  technique.     Let 


«»-eV  *  -t<»-t< 


(14) 


K* 


where 


B*      =       fv,    -    §„     V     &  TV, 

*B»  *  ^  -  P*  ■*  h  ?% 


(15) 


and    |^  ,  VC^  ,  fcC-    are  positive  numbers.    Then  differentiation  of  (14)  with 
times  gives 


-i.| 

u^Wt 

+ 

1 
**> 

Hence, 

,  choosing 

8     4J&+JL.431 


(16) 
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(17) 


results  in  the  equation 


it  3  ~?*tc" 


(18) 


Since  £*<£/  choosing  £    as 


£    »  £  (19) 


or 


S«jvi  £ 


(2  0) 


insures  that  -77  is  negative.     Hence,  we  expect  that  y   vvill  tend  to  zero  and 
the  gains  will  tend  to  the  ideal  values  at  each  flight  condition.    The  system  is 
not  stable  in  any  classical  sense  because  the  manner  in  which  the  gains  con- 
verge depends  on  the  activity  and  frequency  content  of  the  input.     However,  a 
concept  of  eventual  stability  has  been  used  by  LaSalle  and  Rath,  Reference  15, 
to  describe  the  mathematical  behavior  of  such  systems. 

We  prefer  the  mechanization  given  by  (2  0).    Shipley  at  first  used  (19)  and 
found  that  an  automatic  gain-control  circuit  was  needed  to  provide  stability 
under  large  inputs.    The  final  SIDAC  system  (5)  uses  not  only  the  sign  of  the 
error  but  also  the  sign  of  the  dependent  variables  to  form  the  quantities  which 
when  integrated  produce  the  varying  gains.     In  our  system  the  gains  are  found 
from  the  equations 


14 


P,     -     fe,^C^«^ 


(21) 


with  ,  —  S  ._.  —  -.- 

Modifications  to  Include  the 
Dynamics  of  the  Actuator 

Returning  to  Figure  2  we  now  include  the  actuator,  considered  as  a  first 
order  lag,  add  the  corresponding  model  of  the  actuator,  and  the  fixed-feedback 
outer-loop.    The  equations  are  then  modified  to  read 


£        =       G     -   p^©     -  fj„>1   -    p%     T^t  (22) 


Since  the  Laplace  operator  is  equivalent  to  differentiation,  the  first  two  equa- 
tions may  be  rewritten  as 
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Combining  these  two  and  including  the  evaluation  of  the  control  function  (8) 
gives 

(t«s+i)  £  =  TVC^-p^)©    +  t*  ( f ., - p*, ) -^ 

+  i^-H*  &rt)®  *  CP--pH+p*r„)^      <24) 
+  (it  -  n  +  h,r*)c- 

Thus,  the  actuator  introduces  the  first  two  terms  on  the  right  hand  side  which 
were  not  found  in  Equation  12.    These  could  be  taken  into  account  by  adding 
two  more  feedbacks  of   0  and   n    with  varying  gains.    We  avoid  this  compli- 
cation and  hope  that  these  terms  will  not  lead  to  serious  errors.    Analog  simu- 
lation shows  that  this  is  the  case. 

The  requirement  given  by  Equation  5,  that  C*  should  be  a  multiple  of  the 
command  input  is  unrealistic,  considering  the  vast  latitude  by  Figure  1.  We 
modify  this  by  requiring  that  C*  follow  the  command  by  a  simple  lag, 

*                             -It  - 

C         = C-  (25) 

Tc  S  ♦  »  > 

and  calculate  the  deviation  from  this  by 

£ c      =      C         +        Cc    ,  (26) 


This  is  the  equivalent  to 


(27) 
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when  the  definition  of  C*,  Equation  4,  is  used.    The  first  equation  of  (23),  upon 
adding  the  fixed  feedback,  becomes 


(28) 


(29) 


Multiplying  Equation  28  by  a  constant         and  subtracting  it  from  Equation  27 
yields 

We  now  have  the  quantities  Hi   ,   H*,  ,  U«  ,    U    and  K.    at  our  disposal  to  make 
the  C*-error  small.     It  is  assumed  that  the  gain  changers  hold  the  performance 
error  £    small.    Thus,  require 

te  +/*.  ft  -  d 

Uc    +y"P^  «>«ft  *H"° 

The  first  equation  gives  us  a  gain  relation,  the  others  define  the  fixed -feedback 
gains  as 


(3  0) 


U^c^y, 


(31) 
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The  quantity^   can  now  be  taken  as  a  gain  for  a  root  locus  study  or  to  kill  off 
one  of  the  two  remaining  terms  in  Equation  (29).    With  the  data 

Tc     »  .IS-  ?y  =  -3 

the  variations  of  these  gains,  the  open-loop  steady-state  gain  G    ,  and  the  re- 

o 

suiting  closed-loop  roots  are  shown  in  Table  2,   (the  ideal  feedbacks  (13)  are 
assumed).    We  note  that  the  open-loop  gain  is  not  very  large.    Bode  plots  show 
that  the  bandwidth  is  always  less  than  12  radians  per  second.    This  analysis 
was  made  from  the  transfer  functions  recorded  in  the  Appendix. 

Performance  of  the  System 

The  controller  was  evaluated  on  an  all-analog  simulation.    After  initial 
study  with  only  the  actuator  as  the  complication,  a  representation  of  the  pri- 
mary control  system,  a  second-order,  forty-radian-per-second  servo,  a  simula- 
tion of  0.  1-degree  elevator  hysteresis,  and  equations  describing  a  bending 
mode  were  added.    The  system  accepted  these  although  the  rates  and  accura- 
cies for  the  convergence  of  the  gains  were  reduced.    The  dynamics  of  the 
second-order  servo  introduced  limit  cycles  under  relatively  large  command  in- 
puts.    The  effects  of  noisy  measurements  and  air  turbulence  were  studied.    The 
system  was  not  upset  by  these.    The  most  notable  aspect  of  the  performance 
was  the  insensitivity  of  the  C* -response  to  large  deviations  of  the  varying 
gains  from  their  ideal  values. 

The  fixed  outer  feedback  could  be  eliminated  by  choosing  the  parameters 
of  the  error  expression  (11),  |L    ,  0^  ,  (k     to  be  equal  to  the  corresponding 


values  in  the  C*-equation  (16).    The  simulation  showed  that  with  this  choice 
the  gain-changing  mechanism  was  unstable.    To  keep    the  loop-gains  and  the 
errors  introduced  by  neglected  complications  small,  these  parameters  are  picked 
to  be  in  the  range  of  the  corresponding  aircraft  parameters    ft.   ,    By*  /  p*    . 
Then  the  system  is  stable. 

Conclusion 

The  Nonvarying-C*  Control  Scheme  has  been  shown  to  have  practical  appli- 
cation for  fixed-wing  aircraft.     It  should  be  possible  to  extend  the  method  for 
use  on  helicopters.    The  approach  of  choosing  the  variable-gain  mechanism  on 
the  basis  of  system  performance  rather  than  requiring  an  explicit  parameter 
identification  results  in  a  low-bandwidth  system  that  is  very  tolerant  to  noise 
and  high-order  effects.    The  study  also  suggests  a  modification  of  the  C*- 
Criterion  for  longitudinal  handling  qualities. 
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APPENDIX 
Transfer  Functions 


Assume  the  ideal  gains  (13)  obtain.   Then  the  open- loop  transfer  function, 
breaking  the  loop  before  the  actuator,  has  the  numerator 

with  denominator  of  the  aircraft  and  the  actuator  as 
the  open- loop  steady- state  gain  is 

Closing  the  loop  produces  the  denominator 

After  a  little  algebra,  the  C*- transfer  function  is  calculated  to  be 

C*  _   ^<fo   _    (Z*S\   E.StEo  )S 


with 


E,  =    —  \^fy-fyV*r«-- ^»  CTc"t"^tT*fJv,)J 
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